A computational model predicts sex-specific responses to calcium channel blockers in mammalian mesenteric vascular smooth muscle.

The function of the smooth muscle cells lining the walls of mammalian systemic arteries and arterioles is to regulate the diameter of the vessels to control blood flow and blood pressure. Here, we describe an in silico model, which we call the 'Hernandez-Hernandez model', of electrical and Ca2+ signaling in arterial myocytes based on new experimental data indicating sex-specific differences in male and female arterial myocytes from murine resistance arteries. The model suggests the fundamental ionic mechanisms underlying membrane potential and intracellular Ca2+ signaling during the development of myogenic tone in arterial blood vessels. Although experimental data suggest that KV1.5 channel currents have similar amplitudes, kinetics, and voltage dependencies in male and female myocytes, simulations suggest that the KV1.5 current is the dominant current regulating membrane potential in male myocytes. In female cells, which have larger KV2.1 channel expression and longer time constants for activation than male myocytes, predictions from simulated female myocytes suggest that KV2.1 plays a primary role in the control of membrane potential. Over the physiological range of membrane potentials, the gating of a small number of voltage-gated K+ channels and L-type Ca2+ channels are predicted to drive sex-specific differences in intracellular Ca2+ and excitability. We also show that in an idealized computational model of a vessel, female arterial smooth muscle exhibits heightened sensitivity to commonly used Ca2+ channel blockers compared to male. In summary, we present a new model framework to investigate the potential sex-specific impact of antihypertensive drugs.

High blood pressure is a major risk factor for heart disease, which is one of the leading causes of death worldwide. While drugs are available to control blood pressure, male and female patients can respond differently to treatment. However, the biological mechanisms behind this sex difference are not fully understood. Blood pressure is controlled by cells lining the artery walls called smooth muscle cells which alter the width of blood vessels. On the surface of smooth muscle cells are potassium and calcium channels which control the cell's electrical activity. When calcium ions enter the cell via calcium channels, this generates an electrical signal that causes the smooth muscle to contract and narrow the blood vessel. Potassium ions then flood out of the cell via potassium channels to dampen the rise in electrical activity, causing the muscle to relax and widen the artery. There are various sub-types of potassium and calcium channels in smooth muscle cells. Here, Hernandez-Hernandez et al. set out to find how these channels differ between male and female mice, and whether these sex differences could alter the response to blood pressure medication. The team developed a computational model of a smooth muscle cell, incorporating data from laboratory experiments measuring differences in cells isolated from the arteries of male and female mice. The model predicted that the sub-types of potassium and calcium channels in smooth muscle cells varied between males and females, and how the channels impacted electrical activity also differed. For instance, the potassium channel Kv2.1 was found to have a greater role in controlling electrical activity in female mice, and this sex difference impacted blood vessel contraction. The model also predicted that female mice were more sensitive than males to calcium channel blockers, a drug commonly prescribed to treat high blood pressure. The findings by Hernandez-Hernandez et al. provide new insights into the biological mechanisms underlying sex differences in response to blood pressure medication. They also demonstrate how computational models can be used to predict the effects of drugs on different individuals. In the future, these predictions may help researchers to identify better, more personalized treatments for blood pressure.

A computational model predicts sex-specific responses to calcium channel blockers in mammalian mesenteric vascular smooth muscle.

The function of the smooth muscle cells lining the walls of mammalian systemic arteries and arterioles is to regulate the diameter of the vessels to control blood flow and blood pressure. Here, we describe an in-silico model, which we call the "Hernandez-Hernandez model", of electrical and Ca2+ signaling in arterial myocytes based on new experimental data indicating sex-specific differences in male and female arterial myocytes from murine resistance arteries. The model suggests the fundamental ionic mechanisms underlying membrane potential and intracellular Ca2+ signaling during the development of myogenic tone in arterial blood vessels. Although experimental data suggest that KV1.5 channel currents have similar amplitudes, kinetics, and voltage dependencies in male and female myocytes, simulations suggest that the KV1.5 current is the dominant current regulating membrane potential in male myocytes. In female cells, which have larger KV2.1 channel expression and longer time constants for activation than male myocytes, predictions from simulated female myocytes suggest that KV2.1 plays a primary role in the control of membrane potential. Over the physiological range of membrane potentials, the gating of a small number of voltage-gated K+ channels and L-type Ca2+ channels are predicted to drive sex-specific differences in intracellular Ca2+ and excitability. We also show that in an idealized computational model of a vessel, female arterial smooth muscle exhibits heightened sensitivity to commonly used Ca2+ channel blockers compared to male. In summary, we present a new model framework to investigate the potential sex-specific impact of anti-hypertensive drugs.

The formation of KV2.1 macro-clusters is required for sex-specific differences in L-type CaV1.2 clustering and function in arterial myocytes.

In arterial myocytes, the canonical function of voltage-gated CaV1.2 and KV2.1 channels is to induce myocyte contraction and relaxation through their responses to membrane depolarization, respectively. Paradoxically, KV2.1 also plays a sex-specific role by promoting the clustering and activity of CaV1.2 channels. However, the impact of KV2.1 protein organization on CaV1.2 function remains poorly understood. We discovered that KV2.1 forms micro-clusters, which can transform into large macro-clusters when a critical clustering site (S590) in the channel is phosphorylated in arterial myocytes. Notably, female myocytes exhibit greater phosphorylation of S590, and macro-cluster formation compared to males. Contrary to current models, the activity of KV2.1 channels seems unrelated to density or macro-clustering in arterial myocytes. Disrupting the KV2.1 clustering site (KV2.1S590A) eliminated KV2.1 macro-clustering and sex-specific differences in CaV1.2 cluster size and activity. We propose that the degree of KV2.1 clustering tunes CaV1.2 channel function in a sex-specific manner in arterial myocytes.

The formation of K V 2.1 macro-clusters is required for sex-specific differences in L-type Ca V 1.2 clustering and function in arterial myocytes.

In arterial myocytes, the canonical function of voltage-gated Ca V 1.2 and K V 2.1 channels is to induce myocyte contraction and relaxation through their responses to membrane depolarization, respectively. Paradoxically, K V 2.1 also plays a sex-specific role by promoting the clustering and activity of Ca V 1.2 channels. However, the impact of K V 2.1 protein organization on Ca V 1.2 function remains poorly understood. We discovered that K V 2.1 forms micro-clusters, which can transform into large macro-clusters when a critical clustering site (S590) in the channel is phosphorylated in arterial myocytes. Notably, female myocytes exhibit greater phosphorylation of S590, and macro-cluster formation compared to males. Contrary to current models, the activity of K V 2.1 channels seems unrelated to density or macro-clustering in arterial myocytes. Disrupting the K V 2.1 clustering site (K V 2.1 S590A ) eliminated K V 2.1 macro-clustering and sex-specific differences in Ca V 1.2 cluster size and activity. We propose that the degree of K V 2.1 clustering tunes Ca V 1.2 channel function in a sex-specific manner in arterial myocytes.

The formation of KV2.1 macro-clusters is required for sex-specific differences in L-type CaV1.2 clustering and function in arterial myocytes.

In arterial myocytes, the canonical function of voltage-gated CaV1.2 and KV2.1 channels is to induce myocyte contraction and relaxation through their responses to membrane depolarization, respectively. Paradoxically, KV2.1 also plays a sex-specific role by promoting the clustering and activity of CaV1.2 channels. However, the impact of KV2.1 protein organization on CaV1.2 function remains poorly understood. We discovered that KV2.1 forms micro-clusters, which can transform into large macro-clusters when a critical clustering site (S590) in the channel is phosphorylated in arterial myocytes. Notably, female myocytes exhibit greater phosphorylation of S590, and macro-cluster formation compared to males. Contrary to current models, the activity of KV2.1 channels seems unrelated to density or macro-clustering in arterial myocytes. Disrupting the KV2.1 clustering site (KV2.1S590A) eliminated KV2.1 macro-clustering and sex-specific differences in CaV1.2 cluster size and activity. We propose that the degree of KV2.1 clustering tunes CaV1.2 channel function in a sex-specific manner in arterial myocytes.

A stochastic model of ion channel cluster formation in the plasma membrane.

Ion channels are often found arranged into dense clusters in the plasma membranes of excitable cells, but the mechanisms underlying the formation and maintenance of these functional aggregates are unknown. Here, we tested the hypothesis that channel clustering is the consequence of a stochastic self-assembly process and propose a model by which channel clusters are formed and regulated in size. Our hypothesis is based on statistical analyses of the size distributions of the channel clusters we measured in neurons, ventricular myocytes, arterial smooth muscle, and heterologous cells, which in all cases were described by exponential functions, indicative of a Poisson process (i.e., clusters form in a continuous, independent, and memory-less fashion). We were able to reproduce the observed cluster distributions of five different types of channels in the membrane of excitable and tsA-201 cells in simulations using a computer model in which channels are "delivered" to the membrane at randomly assigned locations. The model's three parameters represent channel cluster nucleation, growth, and removal probabilities, the values of which were estimated based on our experimental measurements. We also determined the time course of cluster formation and membrane dwell time for CaV1.2 and TRPV4 channels expressed in tsA-201 cells to constrain our model. In addition, we elaborated a more complex version of our model that incorporated a self-regulating feedback mechanism to shape channel cluster formation. The strong inference we make from our results is that CaV1.2, CaV1.3, BK, and TRPV4 proteins are all randomly inserted into the plasma membranes of excitable cells and that they form homogeneous clusters that increase in size until they reach a steady state. Further, it appears likely that cluster size for a diverse set of membrane-bound proteins and a wide range of cell types is regulated by a common feedback mechanism.

Nonlinear signaling on biological networks: The role of stochasticity and spectral clustering.

Signal transduction within biological cells is governed by networks of interacting proteins. Communication between these proteins is mediated by signaling molecules which bind to receptors and induce stochastic transitions between different conformational states. Signaling is typically a cooperative process which requires the occurrence of multiple binding events so that reaction rates have a nonlinear dependence on the amount of signaling molecule. It is this nonlinearity that endows biological signaling networks with robust switchlike properties which are critical to their biological function. In this study we investigate how the properties of these signaling systems depend on the network architecture. Our main result is that these nonlinear networks exhibit bistability where the network activity can switch between states that correspond to a low and high activity level. We show that this bistable regime emerges at a critical coupling strength that is determined by the spectral structure of the network. In particular, the set of nodes that correspond to large components of the leading eigenvector of the adjacency matrix determines the onset of bistability. Above this transition the eigenvectors of the adjacency matrix determine a hierarchy of clusters, defined by its spectral properties, which are activated sequentially with increasing network activity. We argue further that the onset of bistability occurs either continuously or discontinuously depending upon whether the leading eigenvector is localized or delocalized. Finally, we show that at low network coupling stochastic transitions to the active branch are also driven by the set of nodes that contribute more strongly to the leading eigenvector. However, at high coupling, transitions are insensitive to network structure since the network can be activated by stochastic transitions of a few nodes. Thus this work identifies important features of biological signaling networks that may underlie their biological function.

Role of connectivity and fluctuations in the nucleation of calcium waves in cardiac cells.

Spontaneous calcium release (SCR) occurs when ion channel fluctuations lead to the nucleation of calcium waves in cardiac cells. This phenomenon is important since it has been implicated as a cause of various cardiac arrhythmias. However, to date, it is not understood what determines the timing and location of spontaneous calcium waves within cells. Here, we analyze a simplified model of SCR in which calcium release is modeled as a stochastic processes on a two-dimensional network of randomly distributed sites. Using this model we identify the essential parameters describing the system and compute the phase diagram. In particular, we identify a critical line which separates pinned and propagating fronts, and show that above this line wave nucleation is governed by fluctuations and the spatial connectivity of calcium release units. Using a mean-field analysis we show that the sites of wave nucleation are predicted by localized eigenvectors of a matrix representing the network connectivity of release sites. This result provides insight on the interplay between connectivity and fluctuations in the genesis of SCR in cardiac myocytes.

Publisher's Note: Role of connectivity and fluctuations in the nucleation of calcium waves in cardiac cells [Phys. Rev. E 92, 052715 (2015)].

This corrects the article DOI: 10.1103/PhysRevE.92.052715.